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^ ' Abstract 
Os ; 

. The e~^e~ ujtt^ — > 7r''7r*^7 process was investigated in the SND experiment at 

X ' the VEPP-2M colhder. A narrow energy interval near the i?!)-meson was scanned. 

I ' The observed cross-section reveals, at the level of three standard deviation, the 
interference effect caused by ^ 7r'^7r''7 decay. The cross-section parameters, as well 

r~| , as the real and imaginary parts of the <j)-meson related amplitude, were measured. 
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1 Introduction 

There are several reasons P, [2] , which make the experimental study of the 
e+e" ujTT^ reaction interesting. First of all, it is expected that the radial 
excitations of p-meson should reveal themselves in this process, thus giving 
us a possibility to study their properties. The e+e^ cott^ transition is also 
important for the precise determination of the total cross section of e^e~ 
annihilation into hadrons. In the vector meson dominance model e^e^ 
ujTT^ transition is connected to the topTT vertex, which appears also in a 
number of hadron decays, like u 37r, to 7r^7, p ttj, cu p^ p^n^^ 
TT*^ 27. The precise experimental data about these processes stimulate 
theoretical study of the underlying hadron dynamics. In the low energy 
limit, the chiral perturbation theory [3] and effective chiral Lagrangians [3] 
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give a phenomenological description of meson physics, as it is commonly 
believed nowadays. But their applicability is in question for energies above 

1 GeV. On the other hand, for energies 1-^2 GeV, the perturbative QCD is 
also not applicable. So any precise experimental information in this energy 
range can be considered as " data in searching of the theory" . Simple but 
successful vector meson dominance picture may become insufficient when 
the precision of the data increases. 

The e^e^ utt^ transition can be studied in either e^e^ utt*^ 
Tr^TT^TT^TT^ Or 6^6^ OJTT^ TT^TT^'j chanucls. The former can provide 
about one order of magnitude more statistics, but the latter is more prefer- 
able concerning background conditions. 

The e+e~ lutt^ 7r*^7r*^7 reaction was studied earlier in the en- 
ergy range 1.0 -!- 1.4 GeV by the ND detector P, [2]. Indirectly, the 
cr(e+e" cjTT^) cross section was extracted also from the ARGUS data 
on UrCOTT' decay |[5] , under assumption of the Conserved Vector Cur- 

rent (CVC). The results are in good agreement, thus confirming the CVC 
hypothesis [Bj. 

In 1995 a new set of experiments began with the SND detector |7l El E] 
on the Novosibirsk VEPP-2M storage ring p!0|. The maximum luminosity 
of VEPP-2M at 2E = IGeV equals 3 • 10^^cm-^sec-\ 

Below we report the results on e+e^ utt'^ 7r'^7r'^7 process, based on 
1996-1997 two-year SND statistics. In 1996 several scans were made of the 
energy interval 2E = 984-^ 1040Mey |II]. From the 1997 experiment [12j 
a part of statistics, collected at the center of mass energies 980, 1040 and 
1060 MeV, was used. 

2 Theoretical Model 

As we have already mentioned in the introduction, currently no definite 
predictions can be made from the basic first-principle theory of strong in- 
teractions (QCD) in this energy range. The limited accuracy of the existing 
experimental data leaves enough room for various phenomenologically in- 
spired models. Now a light is seen at the end of the tunnel: several meson 
factories come into operation, and also a new generation of fixed-target 
experiments will accumulate a huge number of events with strongly inter- 
acting particles. So the experimental information is expected to become 
increasingly precise. Of course a common wisdom says that "If we see a 
light at the end of the tunnel, it's the light of an oncoming train" [13]. So 
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it is not excluded that the present theoretical models will be shattered by 
this "train" , but it can deliver a new passenger also. 

Here we consider a simple phenomenological model, based on the vector 
meson dominance scheme. Our aim is twofold: to estimate the expected 
cross section and to establish a framework for Monte-Carlo simulation. 

2.1 General considerations 

Let be the matrix element of the electromagnetic current between the 
vacuum and the tt'^tt'^^ final state. Then the amplitude for the e^e~ 
7r^7r^7 transition via one-photon annihilation diagram 




reads (up to irrelevant phase factor) 

M(e+e" 7rV7) = -v{p+)-f^u{p-)J^ , s = {p++P-f. 

s 

The standard considerations lead to the following cross section with unpo- 
larized beams IH] 



TTQ! 



d(7{e^e~ ^ 7rV"7) = — (Jf • Ji + J| • J2) d^, (1) 

s 

where only transverse (with respect to the beam direction) components of 
J contribute and the invariant phase space element is 

,^ dqi dq2 dk 4 

= (2^ (2^(2^^'"^ - ^1 - ^2 - k). 

should satisfy the current conservation condition Qf^J^ = 0. On the 
other hand = T^^yC", being the photon polarization 4-vector. So T^j, 
is a gauge invariant tensor: 

Q^T^''' = k^T'"" = 0. 
3 



There exists a general procedure how to construct such gauge invariant 
tensors which are free from kinematical singularities |[15] . In our case there 
are just three independent tensors 



7^ 

jJLV 

So 



= • Q)9^iv - k^Q^, Q = qi + q2 + k=p^+p+, q = ^{q^ - q^) 

L^} = {k ■ Q)q^^qly - {k ■ q){k^,qy + q^Q^) + {k ■ qfg^^ 
41^ = (k ■ q)[Q'g,, - Q,Q,] + q,[{k ■ Q)Q, - Q%]. 



dimensionless form factors Ai are determined by the concrete dynamical 
model. Summing over the final state photon polarizations (REDUCE pro- 
gram 117] proves to be useful for these calculations) and performing the 
angular integration, we get [16] 



i (jf • Ji + J| • J2) = -\{TrTi, + rrr2,) XI, X2) = 

Cn\Ai\' + C22|^2|' + Cssl^sl' + CuiA.A; + ^^^2) + (2) 
CuiAiAl + ^t^3) + C23(^2^^ + A;A,), 

where (x = ^ = 2 — xi — X2, = ^, r = ^) 

Cii = ^3:^ , Ci3 = - 0:2) , C23 = ^{xi - X2)^ , 

12 2 

C22 = - X2f + -x'^{r^ -l + xf + -{xi - X2f{l - x)(r2 - 1 + x) 



C33 = ^{xi - X2f{l + x) - -x\r^ - 1 + x) (3) 

2 r. 



C12 = o [(^1 "2:2) + X {r - 1 + x) 
o 

So for the total cross-section we obtain (the factor \ accounts for identical 

Xl+ X2+ 

cr(e+e" 7rV7) = ^ J dxi j dx2 F{2 - xi - X2, xi, X2) 



2567r2E2^ 



/rfx/feF(x,2-x-X2,X2). (4) 
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The integration limits are determined from the condition 

{2E -El- E^f - qi 



cos ^12 = 



'2 -.2 

- q2 



<i, 



which gives 



xi- = r , Xi4_ = 1 , = , x+ = 1 — r' 



X2± = 



1 



1 - Xi + 



-(1 - xijVa^i 



2 y-i2 





\ — r'^ — X 


-l±^ 


1 — X 





(5) 



2.2 Vector mesons contribution 

To proceed, we need expressions for the Ai form factors. It is expected 
that the main contribution comes from diagrams of the following type 




Using kinematical structures of the Vi V27t^ and V vertexes, 
which are determined from the Lorentz covariance, we get for the form 
factors 



A = - 1 + ^xj [g{xi) + g{x2)] + ^{xi - X2)[g{xi) - g{x2)], 

A2 = ^[g{xi) + g{x2)] , As = -^[g{xi) - g{x2)] , (6) 
where 

g{x) = E ^(^^'^^H^) , 
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and individual contributions look like (note that V ^ ttj coupling constant 
is defined as egy-K-f) 

s 1 



gv^ s - M^^ + iMv^Tv, 1 - x + ' 

71/ = -(m^ - + iMvTv). (7) 

To take into account the resonance width dependence on energy, one should 
replace MyTy — ^ y/q^Tv{q^), q being the resonance 4-momentum. 

2.3 Coupling constants 

gv-Kj can be determined from the V 7r^7 decay width. Namely 

^^^^ ~ (M2-mJ)3 • 

Using r(a; — ^ tt^t), r(/9''' — ^ 7r"'"7), and r(0 — ^ 7r°7) as inputs, we get 

g^^^ = (2.32 ± 0.06)GeV-i , gp^^ = (0.73 ± 0.04)GeV-^ , 

g^^^ = (0.138 ± 0.007)GeV-^ . (8) 

As for the photon-vector meson coupling constant, it can be determined 
from the r{y e^e~) width 

g^ ^ o? My 

47T 3 r{V e+e-) ' 

and we get 

gp = 5.0 ± 0.1 , g^ = 17.0 ± 0.3 , g^ = 12.9 ± 0.2 . (9) 

More complex is the situation with g^jpj: coupling constant. Assuming 
vector meson dominance, upir vertex appears in a number of processes. 
For example, V ir^j decay proceeds via 

7 
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and therefore gvn-f = ^^7^- So gwin can be determined from TiV ^ 7r^7) 
and TiVi e^e^). On the other hand, tt^ 2j decay, in the framework 
of the vector meson dominance, goes through 




+ {ki < — > k2) , 



k2 



and therefore 



r(7r0 ^ 27) = Tra V 

\9p9ujj 



So QujpTT can be extracted from r(7r^ 27), T{p e+e ), and r(a; 
e+e-). 

AU these methods for g^^pT^ determination give consistent values: 

(11.7± 0.5)GeV~^ from T{p e+e") and T{uj 7r°7) , 
(12.6 ± 0.7)GeV"^ from r{u e+e") and r{p+ 7r+7) , (10) 

(12.5 ± 0.9)GeV"^ from r(7r° ^ 27), r{uj 7r%), r{p+ tt+'j) , 

(11.8 ± 0.5)GeV"^ from r(7r° ^ 27), r{io e+e"), r{p e+e") . 

But if the value of gi^^pj^ is extracted from the experimental to ^ Stt decay 
width, assuming that this decay proceeds through to ^ pir intermediate 
state and that the gpT^-j^ coupling constant is determined from the r{p ^ 
27r), one obtains 

= (14.3 ± 0.2) GeV-^ (11) 
This is closer to the chiral model prediction [18] 



9upn = « 14.9 GeV^i , A « 93MeV 



gcup-K can be estimated also from QCD sum rules |I9], [20] with the result 
g^p^ = (16 ±2)GeY-' [20]. 



In Monte-Carlo simulation, we have used ( fTTI ) and (/^pTr ~ — O.SlGeV"^, 
which follows from r(0 Stt). These values assume only p-pole decay 
mechanism for cu ^ Sir transition, without a possible cj ^ Stt contact 
term [211 • iact, the uncertainty in these coupling constants does not 



effect significantly the Monte-Carlo estimates for the detection efficiencies. 
2.4 (f)—p—u mixing contribution 

(J) UJ7T vertex is forbidden by G-parity, which is negative for all three 
particles. Another way to see this is to notice that it is impossible to 
construct isospin invariant trilinear coupling between and cu isosinglets 
and isovector pion. 

But the isospin symmetry breaks due to electromagnetic effects and 
mass difference between u and d quarks. As a result, pure isospin eigen- 
states uj and pi mix and cease to be mass eigenstates. Instead we will 
have a nondiagonal mass matrix 



( Pi 



4 z^p \( Pi 



Mass eigenstates (physical p and uj mesons) are mixtures of isospin eigen- 
states, which for small p-u mixing angle e^jp ( sin e^^p « e^^p) look like 

p = Pi - e^pUi 
u = ui + e^ppi 

The mixing angle e^^ is determined from the condition 

Zp Zojp\ ( Pi \ ( \ ( Zp \ ( p 



CO 



Zu;p Z^ I \ Ui I ^ M / \ CJ 



which gives in the first order in z^jp and e^jp [22 



P ^ '^p 1 '^UJ ^ '^UJ 1 ^LUp 



Zu) Zp 



Here [221 zy = {My — ^^)^ ~ My — iMyTy is the resonance complex mass 



square. So 

e = ^ (12' 



As for p-LU mixing amplitude z^p, we will take the value 

z^p = (-3800 ± 370)MeV'^ , (13) 
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which follows from the pion form factor studies [2^. 

Since pi = p + t^pOJ^ the presence of the piir vertex will induce an 
effective ujtt vertex with a coupling constant g^^J!^ « ^ujpg(f>p'n ■ 

The described above picture of p-cu mixing can be trivially generalized 
to include mixings with the (f) meson also [23] : 



P = Pi — ^.yjpLOi — €(f>p(pi 
cj = cj/ + Cujppi - e^(f,(j)i . 
(J) = (f)i + e^ppi + euj4)Ui 

Due to (j)-p mixing another contribution to the (j) cutt transition will arise 
9^(^Sn ~ ^<l)p9ujpTT • Because Igcupnl ^ bi^ipTrl, this latter contribution turns out 
to be of the same order as g^^j}. So for the ^f^^^Tr coupling constant we will 
use 

9(f>U>Tr 9(f)U!TT ~l~ 9(I)UJTT ^bJp94'pTt ~^ ^4}p9ujp-K • (-'-^) 

The corresponding contribution in e+e^ 7r*^7r^7 is then determined by 
equation (7). 

The 0-/9 mixing parameter e^p can be extracted from the e^e^ ti^t(~ 
experimental data near the 0-meson j23j : 



- z (15) 

where Z is experimentally measured interference magnitude pi|. Taking 
average values Re{Z) = (8 ± 2) • IQ-^, Im{Z) = -(3.5 ± 1.3) • 10"^ from 
15, we obtain 

e^p = (0.72 ± 0.18) • 10"^ - z(0.87 ± 0.32) • 10~^ . (16) 
2.5 p' meson contribution 

No definite theoretical predictions exist for /9'(1465) meson contribution, 
nor it is definitely known whether the only one radial excitation gives 
significant contribution in this energy region. We will use the following 
parameterization for this part of the g{x) function 

g^P''^\x) = Re'^k ^ ^ 



^s-M'i + iMp^Tp' l-x + 7^ 



lu. = -{m^-M^ + iM^r^), (17) 
s 
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where kp = iraM^ ^^"''^^^^^ corresponds to the /o-meson contribution. R and 

^ parameters were determined by fitting a{e^e^ utt) = ^^^^(^!^o^^ 
cross section to higher energy experimental data from . The fit gives the 
following values: 

i? = 0.72 ±0.03, <e = 3.25 ± 0.08 . (18) 

Our results of this fit confirms the conclusions of [25] that the /9'(1465) 
meson contribution is sufficient to describe the existing experimental data 
in the 1.05 1.6 GeV energy range. 

Note that we use the PDG |233 values Mp> = 1.465 GeV and T^/ = 
0.31 GeV. To simulate roughly a threshold effect due to the dominant ptttt 
decay of the p' meson, we have assumed, as in [2S], that the p'- meson width 



rises linearly from = 0.8 GeV up to = Mp>, remaining constant 
above. 



2.6 Theoretical prediction 

Using the estimates for various coupling constants given above, the total 
cross section can be evaluated for e+e" reaction. The result 

is shown on FiglH Note the significant interference effect near meson, 
although the 0-meson mediated amplitude itself remains small and corre- 
sponds to i?r(0 piT ^ 7r^7r^7) « 10^^, in consistence with [271. The 
(p-p-cj interference contribution into cjtt transition corresponds to 
Br(0 utt) ~ 5.7- 10~^ in excellent agreement with the recent experimen- 
tal result [28] (4.8i j;7 ± 0.8) • 10"^, although it should be beared in mind 
that there may exist also other sources for this transition [29] 



3 Detector and experiment 

SND is a general purpose nonmagnetic detector. This new detector com- 
bines advantages of its predecessor ND detector [3D] and famous Crystal 



Ball detector [31] , that is a good uniformity over the solid angle due to the 
spherical shape, a good e/ir and separation due to multilayer struc- 

ture of the electromagnetic calorimeter, and high hermeticity. The main 
part of the SND detector (Figj2]) is a 3-layer, spherical, highly granulated, 
Nal(Tl) electromagnetic calorimeter [32]. Tracking system, located in the 



detector center, consists of two coaxial cylindrical drift chambers. In the 
radial gap between them a 5 mm thick plastic scintillation counter with a 
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wavelength shifter fiber readout is placed. From outside the calorimeter 
is covered by a thick (12 cm) iron absorber, which attenuates remnants of 
electromagnetic showers. An outer muon/veto system is located outside 
the absorber. It consists of sheets of plastic scintillator and streamer tubes. 

The detailed description of the SND detector can be found elsewhere 
[71, El E]- Here we repeat only the main points, relevant to this study. 

3.1 Calorimeter 

Each layer of the calorimeter includes 520-560 crystals of eight different 
shapes. Most of the crystals have shapes of truncated tetrahedral pyramids. 
The total solid angle covered by the calorimeter is equal to 0.9 of 47r. The 
remaining space is occupied by magnetic structure elements of the storage 
ring, mainly by quadrupole lenses. Pairs of counters of the first two layers 
with thickness of 3Xq and 5Xq respectively are sealed together in common 
containers made of thin (0.1 mm) aluminum foil. In order to improve light 
collection efficiency and to separate one crystal from another, each crystal 
is wrapped in aluminized mylar. The gaps between adjacent crystals do 
not exceed 0.5 mm. The r.m.s. value of the nonuniformity of the light 
collection efficiency along the crystals is less than 3 %. All the containers 
are fastened to 5 mm aluminum supporting hemispheres. The outer layer 
of 6Xq thick counters has a similar design. Thus the total thickness of the 
calorimeter is 13.5Xo (35 cm) of Nal(Tl). 

As photosensitive devices for the calorimeter counters vacuum phototri- 
odes are used [33j . The quantum efficiency of their photocathodes is about 
15 %, average gain is 10, and light collection efficiency for individual coun- 
ters is about 10 %. Signals from phototriodes are amplified by charge 
sensitive preamplifiers located directly on the counters. Output signals are 
carried to shaping amplifiers via 20 m long twisted cables. For the trig- 
ger needs the calorimeter crystals are logically organized into "towers" . A 
tower consists of counters located within the same 18° interval in polar and 
azimuthal directions in all three layers. The number of counters in a tower 
is 12 at large angles and 6 in the regions close to the beam. In addition 
to signals from individual counters each tower produces an analog total 
energy deposition signal and two trigger signals. In order to equalize con- 
tributions from different counters into the total energy deposition signal 
and to obtain equal energy thresholds for trigger signals over the whole 
calorimeter, all shaping amplifiers are equipped with computer controlled 
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attenuators, allowing to adjust channels gain in steps of 1/255. The result- 
ing electronics noise is close to 0.3 MeV (r.m.s.). The dependence of the 
calorimeter energy resolution on photon energies was fitted as [31] 



(Je/E{%) = , 7 ^ . (19) 



The calorimeter energy resolution is determined mostly by the amount 
of passive material inside the calorimeter. Noticeable contribution comes 
also from the shower energy leakage due to limited thickness of Nal(Tl) 
and gaps between crystals. It turned out that the nonuniformity of the 
light collection efficiency in the Nal(Tl) crystals of the inner calorimeter 
layer contributes significantly into energy resolution. The simulation of 
calorimeter response agrees reasonably well with experiment after these 
effects are included into simulation program. 

Absolute energy calibration of the calorimeter is performed by using cos- 
[35] and Bhabha scattering j31|. The latter process, together 



mic muons 

with the two-photon annihilation reaction, was also used for the luminosity 
measurements. The systematic error in the integrated luminosity determi- 
nation is about 3%. 

After the calorimeter calibration with e^e^ e^e^ events, the photon 
energies turned out to be biased by about 1%. In order to compensate for 
this bias, the calibration coefficients for photons were corrected accordingly. 



3.2 Tracking system 

The tracking system consists of two cylindrical drift chambers. The length 
of the chamber closest to the beam is 40 cm, its inner and outer diame- 
ters are 4 and 12 cm respectively. The corresponding dimensions of the 
outer chamber are 25, 14 and 24 cm. Both chambers are divided into 20 
jet-type cells in azimuthal plane. Each cell contains 5 sense wires. The 
longitudinal coordinate is measured by the charge division method with 
an accuracy of 3 mm. In addition, a cathode strip readout for inner and 
outer layers provides the improvement of the latter value to 0.5 mm. The 
outer drift chamber improves pattern recognition for multiparticle events. 
The overall angular resolution of the drift chamber system is 0.5 and 1.7 
degrees in azimuthal and polar directions respectively. The impact param- 
eter resolution is 0.5 mm. The solid angle coverage for the inner chamber 
is 96 % of 47r. 
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3.3 Muon detector 



The SND muon system consists of plastic scintillation counters and strea- 
mer tubes [36]. It is intended mainly for cosmic background suppression 
for the events without charged particles. The probability of its triggering 
by the events of the reaction e^e" 77 at the maximum available energy 
is less than 1 %. Muons of the process e^e^ A^-^A^ 5 starting from the 
beam energy about 500 MeV, penetrate through the absorber and hit the 
muon detector. 

3.4 Trigger 

The three-level trigger of the detector selects events of different types: 
events with photons only, events with charged particles, and cosmic muon 
events for the detector calibration. The drift chamber first level trigger 
(FLT) searches for tracks in the drift chambers with impact parameter 
Ar < 20 mm. The calorimeter FLT uses the total energy deposition in the 
calorimeter and double coincidences of calorimeter towers with a threshold 
of 30 MeV. 

Detector electronics provides also signals for a second level trigger. But 
at present the second level trigger proved to be unnecessary and was not 
implemented, although a slot is left for it in the electronics. 

The third level trigger (TLT) is implemented as a special fast computer 
code on a main data acquisition computer. It checks events before recording 
them on tape, rejects cosmic events with charged trigger, and suppresses 
beam-related background using ^-coordinates of tracks, measured by drift 
chambers. In addition TLT identifies coUinear events of Bhabha scattering 
and 27 annihilation. These events are used for monitoring the collider 
luminosity. 

3.5 Data acquisition system 

SND data acquisition system is based on the fast electronics modules 
KLUKVA [37], developed in BINP specially for purpose of the detectors 
CMD-2, KEDR and SND. 

Analog signals from the detector come to the front-end amplifiers and 
shapers, located near the detector. Then the signals are transmitted via 
screened twisted pair cables to the digitizing modules (ADC, TDC, etc) in 
KLUKVA crates. Logical signals from discriminators in KLUKVA crates 
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are collected by the first level trigger interface modules. After 1 /is, needed 
for FLT decision, FLT generates a signal to start digitizing. The contents of 
KLUKVA digitizing modules are read into the RAMs of Output Processor 
modules (OP), located in KLUKVA crates. This procedure finishes in 
120 lis. There are two independent RAMs for events in each OP and one 
event can be stored in digitizing modules, such a fast buffering greatly 
decreases dead time. The total dead time in KLUKVA is 200 /is per event 
selected by FLT. The data from OP RAMs are read by the VAX server 3300 
through the CAM AC interface modules with the rate of 2 ms/event. The 
TLT program can process events with a maximum rate of 45Hz. Finally 
the processed events are put on 8 mm 5 GB EXABYTE tapes. 



3.6 Experiment PHI96 



A description of the PHI96 experiment was published in |[TT|. The PHI96 



experiment was carried out from February 1996 up to January 1997. Seven 
successive scans were performed in 14 energy points in the range 2E from 
980 to 1044 MeV. 

The data sample, which was analyzed, corresponds to an integrated 
luminosity of 4.5 ^6~\ collected by SND in the narrow energy interval near 
the 0-meson. Estimated number of produced 0-mesons equals to 8.3 • 10^. 



4 Event selection 



While studying the channel 

e+e" ^ UTT^ 7r\% (20) 

one should be aware about the possible background from the following 
processes 

e+e" ^ (t)^T]-i ^ 37r% , (21) 

e^e" ^ ^ KsKi neutral particles . (22) 

A primary selection of the e^e^ cjtt*^ 7r^7r^7 candidates was done 
according to following criteria: 

• the event must contain exactly 5 photons in the calorimeter and have 
no charged particles, 

• the azimuthal angle of any final photon lies within the interval < 
< 153^ 
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• the total (normalized over 2E) energy deposition of final photons is 
in the range 0.8 < Etot/^E < 1.1, 

• the normalized full momentum of the event {Ptot/'^E) is less than 0.15. 

The two latter conditions eliminate the main part of the background, orig- 
inated from the K^Ks decay, but they do not help much against 
another background, coming from the process (1211). 

It is interesting to check whether we really have two tt^-s in our 5-photon 
events. It is unlikely for the most energetic photon to come from ti^ decay in 
either (12U|). ( 1?T| ) or e+e^ KlKs 27T^Ki reactions. It is also not likely 
for TT*^ to produce the two softest photons. If we discard the correspond- 
ing combinations, only two possibilities (2, 4), (3, 5) and (2, 5), (3, 4) are left 
for the photons to compose tt^'s (photons are arranged according to their 
energy, the most energetic being the first one). A combined 2-dimensional 
plot of invariant masses of these photon pairs (M24 versus M35 plus M25 
versus M34) is shown in Fig|3l We clearly see that our 5-photon events are 
predominantly 27r^ events also. Besides, FigJ2] illustrates that background 
(EID produces a wider distribution. This can be used for the rejection 
of this background. For this purpose, a kinematic fit was performed for 
each 5-photon event under the assumption that there are two tt^'S in the 
final state and energy-momentum balance holds within the experimental 
accuracy. A of this fit (xIot^o^) can be used for the background rejection. 

Background from (EID simulates events mainly due to loss of pho- 

tons through the openings in the calorimeter around detector poles and/or 
merging of close photon pairs. When photons merge in the calorimeter, the 
corresponding electromagnetic shower IS, as a rule, broader in transverse 
direction than the electromagnetic showers from the single photons. This 
circumstance can be used to discriminate against merged photons and so 
against a great deal of background (EI])- The corresponding parameter 
{(j) is described in [SBl- A 2-dimensional distributions of our events in the 



xla^Q^, C7 plane, as well as Monte-Carlo simulated (EDj) signal events and 
(EID background events (FigJlj), indicate that our signal events are almost 
completely bound in the C7 < 20, xIottO-/ < ^0 area. This was confirmed 
by experimental events outside of the 0-meson, where there should be no 
background from the process (I2T|). 

On the basis of these considerations, we have chosen the following two 
sets (Cut I and Cut II) of selection criteria for the channel ( 12U| ) separation 
(in addition to the primary selection rules, described above): 
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the normalized full momentum of the event is less than 0.1 (for Cut 

I). 



• there are two 7r^-mesons in the event, that is, one can find two distinct 
pairs of photons with invariant masses within ±30 MeV from the tt^ 
mass. 



X^Ot^o^, the of the kinematic fit, is less than 20 for Cut I, or is less 
than 40 for looser Cut II, 



• C-^, the parameter describing the transverse profile of the electromag- 
netic shower, is less than for Cut I, or is less than 20 for Cut II. 



One more source of background is a recently observed process |[I2l EH 

e+e" ^ ^ /o7 ^ 7r\% . (23) 

The recoil mass of the photon from this process is peaked at the /o-meson 
mass and this peculiarity can be used to separate a great deal of such 
events from the events of process (1201). We have chosen the My < 700Mey 
condition as one more cut to select events from the process (EUf). where 
stands for the photon recoil mass. 

After applying these cuts, the cj-meson peak is clearly seen in the in- 
variant mass of TT^ and 7 (for each tt^tt^j event, from two possible (7r^,7) 
combinations, the one is taken, which has M^^o^ closest to M^), as it is il- 
lustrated by Figj5] (Cut I). Finally, to extract channel (EOl), the 750MeV < 
Mt^o^ < 820MeV condition was added to the above mentioned cuts. 

The Mj^Qry histograms were fitted by a function 



X-P2) 



2 



Pi exp( r-^ — ) +Pa ■ 

The fit indicates a good agreement between the experiment and simulation. 
For the 2356 MC events , which passed Cut I, the fitted value of the 
constant term is 2?4 = 9.2 ± 1.3. Since in the 750MeV < M^o^ < 820MeV 
interval we have 35 histogram bins, such background constant corresponds 
to 322 ±46 combinatorial background events, that is about 14%. The same 
conclusion follows also for Cut II. 

For the 560 experimental events which passed Cut I, the following back- 
ground constant was obtained: ^4 = 3.3±0.9, which correspond to 116±32 
background events. Subtracting the expected combinatorial background of 
76 ± 11, we can estimate the residual background from (1211), (I22l), (I23|). as 
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40 ± 34 events (~ 7%). For Cut II this background was found to be twice 
higher, that is 133 ± 38 background events from the selected 864 experi- 
mental events (~ 15%). 

On the other hand, Monte Carlo studies indicate the following proba- 
bilities for events from the main background sources to pass Cut I: (2.7 ± 
0.3) • 10-4 for process (EH), and (2.9 ± 0.2) • 10"^ for process (ESI). This 
corresponds to 28 ± 3 background events from process ( 1?T| ) and 22 it 4 from 
process (I23|). So the total number of background events from these sources 
is estimated to be 50 dz 5. For Cut II we expect 116 ± 8 and 27 ± 5 back- 
ground events, respectively, or in total 143 ± 10 events. These numbers are 
in a good agreement with those obtained from the above mentioned fitting 
of the a;-meson peak. 

To check the accuracy of Monte Carlo simulation the following proce- 
dure was applied. For the above described Cut I and Cut II, one of the 
parameters is released and the resulting experimental distribution Hexp 
for this parameter is compared with the sum of MC predictions for the 
signal H^T^f) and expected backgrounds from ( j?Tl ) and ( 1^ sources: Hmc = 
kuj-K^H^T^o + krj^Hjjj + kfo^Hfo^. The normalization coefficients krj-y and kpry 
are determined by the total statistics of generated MC samples for this 
reactions. The normalization coefficient for the signal k^^T^o is then deter- 
mined from the condition that Hexp and Hmc histograms have the equal 
total numbers of events. Note that any variation in the k^-j^o coefficient in- 
dicates some systematics and/or other background sources, not accounted 
for in the comparison. For Cut I the averaged deviation turned out to be 
3% and for Cut II - 4%. 

The distributions for all parameters used in event selections show good 
agreement between MC and experiment. This indirectly indicates that the 
possible background from process (1221 ) is rather small and does not exceed 
~ 5%. As an example, on Fig|S]we present M^^o^ distributions. 

5 Data analysis and results 

We assume the following parameterization for the visible (detection) cross 
section cr„: 

a, = e[l + 5{s)]a{s) + kbaeis) , (24) 

where e is the detection efficiency for the process (I20|). 6{s) accounts for 
the radiative corrections, which are calculated according to the standard 
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procedure [40], (Jb{s) is a background cross section, which was assumed 
to coincide with a{e^e~ ^ 777 ^ 37r%), k is the background suppression 
factor, h is described below. Finally, cr(s) is a cross section of the process 
under investigation. Because we are interested in cr(s) in a narrow energy 
interval, we have taken 

<j{s) = [ao + a'i2E - M^)] \R\^ R=l- Z , ^ . (25) 

The detection efficiency e was calculated using Monte-Carlo simulation in 
conditions of individual scans for various energies. Detection efficiencies do 
not show any significant energy dependence. So we have taken an efficiency 
value averaged over scans and energies, e = (29. 7± 0.25)% (only statistical 
error is indicated) fair estimate for Cut II. 

For a tighter Cut I some systematic errors could be expected. To esti- 
mate this systematics, we compared the numbers of rejected events for each 
parameter of Cut I in the above described Hexp and Hmc distributions. 
The following correction factors {eMc/<^EXp) were obtained: 1.032 ± 0.018 
for Ptot/^E, 1.015 ± 0.023 for xjo^o^, 1-074 ± 0.028 for C7, 1-015 ± 0.023 
for and 0.98 ± 0.028 for M^^o^. In total, the Monte Carlo simulation 
overestimates the detection efficiency 1.12 ± 0.06 times, if we assume that 
there are no correlations between the used selection parameters. With 
this correction factor taken into account, an average detection efficiency 
e = (20.7 ± 1.1)% was obtained for Cut I. 

The background suppression factor k was also assumed to be energy 
independent. It was calculated using ~ 2.15 • 10^ simulated events from 
the process (ED and equals (2.7±0.3) • 10"^ for Cut I and (1.10±0.06) - lO^^ 
for Cut II. 

As we have seen above, for Cut I a relatively large fraction of background 
is expected from process (ESj). Some small amount of ^ KsKi decay 
related background ([22]) is also not excluded. To take into account these 
and other 0-meson related backgrounds, the factor b is introduced in (I24l). 
It is assumed that the different energy dependencies of various background 
cross-sections is not relevant at the present level of statistical accuracy and 
so they all can be approximated by the cr(e+e^ 777) behavior. 

The fit results are given in the Table [B 

Because of specific behavior of the interference effect, it can happen that 
the resonant background cancels the interference dip and so mimics the no- 
interference no-background situation, making it impossible to distinguish 
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parameter 


Cut I 


Cut II 


ctq (nb) 
a' (nb/MeV) 
Re(Z) 
\m{Z) 
h 


0.61 ± 0.06 
(0.43 ±0.24) • 10-2 
0.1 ±0.1 
-0.19 ±0.15 
2.8 ± 1.9 


0.63 ± 0.05 
(0.49 ±0.21) ■ 10-2 
0.12 ±0.08 
0.05 ±0.08 
0.1 ±0.3 


xVd.f. 


11.0/10 


14.5/10 



Table 1: Fitted parameters for the fit with the interference. 

between these two options during the fit. This is confirmed by the following 
observation: If we assume no interference effects in the a{s) and just take 
a{s) = do + a'{2E — Mff^), this linear fit will result in the Table [21 



parameter 


Cut I 


Cut II 


ctq (nb) 
a' (nb/MeV) 
b 


0.58 ± 0.05 
(0.44 ±0.23) ■ 10-2 
0.2 ± 1.1 


0.64 ± 0.03 
(0.58 ±0.20) ■ 10-2 
0.2 ±0.2 


xVd.f. 


12.1/12 


16.7/12 



Table 2: Fitted parameters for the linear fit. 

In spite of our considerations above, from which one could expect 6 ~ 1 
for Cut II (and 6 « 2 for Cut I - the expected background is indeed domi- 
nated by the process ( 1^ for Cut II, while for Cut I half of the background 
originates due to process (I23l)). this fit indicates no-background for Cut II, 
so in fact providing an indirect evidence in favour of the interference effect! 

As we see, the results for Cut I and Cut II are consistent after the 
detection efficiency for Cut I is corrected against estimated systematics. 
Of course, we do not know if there is some part of systematic errors left 
uncorrected for Cut I. On the other hand, for Cut II more background is 
expected and we can neither very precisely estimate this background (for 
example, the part coming from (1221)) nor subtract it during fit (note that 
all of the above fits gave very large errors for the background constant b). 

Therefore we select as a fair estimates for the ctq and a' parameters the 
mean values between Cut I and Cut II, and as systematic errors we take 
the difference between them. Another part (~ 3%) of the systematic error 
can arise from the luminosity measurement. 

One more source of systematics is a possibility that a sixth false photon 
is piled up on the five photon event due to the beam related background. 
This effect is hard to simulate, so the MC efficiency is expected to be 
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somewhat overestimated. To find the correction factor, the 6-photon events 
were investigated as described in the next section. This factor turned out 
to be 1.04 ± 0.02. Our final results in this analysis, which include this 
correction factor, as well as the estimated total systematic error, are 

do = (0.64 ± 0.08) nb 

a = (0.48 ± 0.27) • 10"^ nb/MeV (26) 

As for Re{Z) and Im{z)^ the still high level of background precludes an 
estimation of these parameters. For Cut I, which is more pure against 
background, Im{Z) = —0.19 ± 0.15 is only one sigma effect. 

We have also tried a{s) = C(Jth{s) in a fit, where (Jth{s) is the theoreti- 
cal prediction (including /o'-meson) for the reaction (1^. discussed at the 
beginning. The fitted parameters are given in the Table 3. 



parameter 


Cut I 


Cut II 


c 

h 


1.19 ±0.09 (1.54 ±0.12) 
2.7 ±1.0 (2.8 ±1.0) 


1.29 ±0.08 (1.68 ±0.11) 
1.1 ±0.3 (1.2 ±0.3) 


xVd.f. 


11.7/13 (11.7/13) 


17.5/13 (17.4/13) 



Table 3: Fitted parameters for the fit with the theoretical cross-section. The parameters 
for the fit without p' meson are also given in the parenthesis. The indicated ^-values 
should be multiplied by the correction factor 1.04 ± 0.02, as discussed in the text. 

The obtained cross section is close enough to the theoretical prediction, 
when p'-meson contribution is also included, although the experimental 
cross section is somewhat (~ 25%) higher. This can indicate that more 
then one radial excitation of the p-meson contributes or that the used 
higher energy ee utt cross section is subject to some systematics, which 
we do not take into account while extracting the p' related parameters R 
and ^. Finally, our parameterization for the p' contribution can also be 
not quite adequate (for example, we had only roughly modeled the energy 
dependence of the p' width). In any case, the experimental cross section 
can not be explained by only the p-meson tail, which gives about 1.6 times 
lower result. 

6 Observation of the interference effect 

During the study of the reaction e+e^ utt^ tt^tt^tt'^tt^ , the (j) utt^ 
decay was observed for the first time with the branching ratio about 5 • 10^^ 
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2B]- The decay reveals itself as an interference wave on the nonresonant 
cross section of the process e^e~ cott^. In principle, a similar picture 
should be observed in neutral channel (EDI), as was explained above in 
the theoretical introduction. Really the whole situation here looks more 
complicated because of other meson neutral decays like (J) p^ir^, (J) 
/o7, (77 [32], which have the same final state and interfere with the process 



The interference amplitude with the decay is expected to 

be about 10%, which is close to the value 17% due to ^ cott^ decay, 
obtained in [28]. In our preceding study |12] of the reaction (1201). we did 



not observe the interference because of small statistics and nonresonant 
background. 

In the present work the analysis given above also does not reveal the 
interference effect with certainty in spite of higher statistics, because the 
background is still high. Here we present another analysis specially dedi- 
cated to the interference observation. 

The events selected by the criteria, described above, were divided into 
following 3 classes: 

1) xk"7<20, |M,o^ - 782| < 30, C7 < "5; 

2) xko-y < 20, 30 < |M,o^ - 782| < 60, C7 < "5; 

3) xlo^o^ < 20, |M^o^ - 782| < 30, = 6. 

The 6-photon events were put into the last class. It was done to investigate 
a probability that sixth false photon is superimposed on the five photon 
event due to beam-related background. In the kinematic fit one of photons 
was supposed to be spare. 

The simulation shows that the signal to background ratio is maximal 
in class 1, which was used later for investigation of the interference effects. 
The contribution of process ( 1201 ) in class 2 is about 10 times lower, which 
allows to extract from the data the resonant background from the above 
mentioned processes with , KsKl and 777 intermediate states. Using 
the ratio of the resonant background in classes 1 and 2 obtained from the 
simulation, one can estimate the background in class 1, when the resonant 
background in class 2 is determined from the experimental data. The ratio 
of the background event numbers in the regions \M-,^o-y — 782 1 < 30 and 
30 < iMjrO^ — 782 1 < 60 weakly depends on the hmits on x^^-k^-^ and C7 
parameters. For the process ( I23| ) this ratio varies in the limits 1.09 1.12, 
while for the process (EID — from 0.9 to 1.07. These variations are within 
the statistical accuracy of the corresponding MC simulation. The value 
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I. 1 ± 0.2 was taken for this ratio. 

The visible cross section in each class was presented in the following 
form: 

ayis = a^e[l + 6{s)]a{s) + P^asis) 

where a{s) is given by ([25]), £ is the total over all classes detection efficiency 
for the e'^e~ utt^ tt'^tt^^ process, and 6{s) represents the standard 
radiative corrections [|40j. As above, aB{s) is assumed to coincide with 
a{e~^e~ ^ 777 ^ Stt^j). The parameters ai are the probabilities for our 
process to be found in i-th class, the parameters Pi represent nonresonant 
cross section in i-th class, normalized to the e+e~ ^ ^ 777 process 
cross section. As it was already mentioned, the meson excitation curve 
was described by the main background process (1211). The differences in 
excitation curves for different background processes were neglected at the 
present level of accuracy. The fitting of the visible cross section was done 
for all three classes simultaneously. The parameters ctq, ct', Re{Z), Im{Z), 
Pi except Pi and ai were free. Pi were found from the expression 
Pi = (1.1 =b 0.2) • P2, Oil ~ from normalization = 1- The detection 
efficiency £=27.5%, obtained from the simulation, does not depend on the 
energy. The coefficients di are also energy independent. The total number 
of fit parameters is 8. In each class the cross section was measured in 15 
points. The following values of interference parameters were obtained 

Rt{Z) = 0.036 ±0.052, 
/m(Z) = -0.186 ± 0.063. (27) 

The statistical errors for this parameters are much higher than the system- 
atical ones. The quoted errors include the systematics due to background 
subtraction in class 1. 

The ratio 0^3/ (ai + 0^3) gives the probability that the signal event will be 
lost because the additional background photon is superimposed on it. This 
probabihty turned out to be 0.04 ± 0.02 and the corresponding correction 
factor was included in the ctq and o' determination in the previous section. 

The detected cross section for class 1 and the fitted curve with /d.f. = 

II. 9/11 are shown in Fig|7l The fitted resonant background is also shown 
at the bottom. One could see, that in spite of imposed strong cuts in class 
1, the resonant background is about one third of interference amplitude 
wave and is the dominant source of systematic error in Z . This systematic 
error is estimated to be about 6%. 
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7 Conclusions. 



In conclusion, we obtained in the present work the following values of the 
e+e^ UTT^ 7r'^7r^7 process cross section parameters: 

ao = (0.64 ±0.08) nb, 
a' = (0.0048 ± 0.0027) nb/MeV, 
Re{Z) = 0.036 ±0.052, 
Im{Z) = -0.186 ± 0.063. (28) 

The measured interference amplitude is three standard deviation above 
zero. The measured in [28] interference amplitude for the decay ujtt^ 
is Re{Z) = 0.104 ± 0.029, Im{Z) = -0.118 ± 0.031. Our calculation 
of the interference amplitude for the decay pir'^ tt'^tt'^^, in the 
framework of the model described in the theoretical introduction, gives 
Re[Z) = —0.079, Im{Z) = —0.053. The sum of these contributions 
Re{Z) = 0.025, Im{Z) = —0.171 agrees with our measurement. Otherwise 
stated, if one subtract from the interference amplitude, measured in this 
work, the expected contribution from the (f) p7r^ transition given above, 
the (j) L07T decay branching ratio can be estimated to be 5.4 • 10^^ with 
almost 100% errors, which should be compared to (4.8i|;7±0.8) • 10^^ from 
[2B1 • Note that we did not make a theoretical estimate for contribution into 



interference amplitude from the processes (f) /07, (17 —> 71*^71^^. More 
experimental information about these transitions is required to estimate 
such contribution correctly. 

The measured nonresonant cross section at 0-meson cr(e+e~ uir^) = 
cr(e+e~ ujTT^ TT^TT^^)/ Br{u 7r%) = (7.5 ± 0.9) nb agrees with the 
result (8.7 ± 1.0 ± 0.7) nb from H and with the result (8.2 ± 0.2 ± 0.9) nb 
from [2B1 in channel with charged pions e~^e~ utt^ tt^tt^tt^tt'^ , as well 
as with the preliminary CMD-2 result [Sj (7.2 ±0.5) nb. Such nonresonant 



cross section can not be explained by only /o-meson contribution, which is 
about 1.6 times lower. Although the inclusion of the p' {14:65) meson in 
the fit improves somewhat the agreement between the theoretical predic- 
tion and the experiment, the question about p-meson radial excitations 
lies beyond the scope of this work, because it requires an experimental 
information for higher energies. 

At last, in Table 4 we provide numerical values of the measured cross 
section cr(e+e^ lott^ tt^tt^^) = ^j^^^jj- 
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Vs, MeV 


cross section, nb 


984.1 


0.57 ± 0.11 


1005.0 


0.56 ± 0.14 


1010.7 


0.47 ± 0.10 


1016.2 


0.55 ± 0.10 


1017.1 


0.48 ± 0.10 


1018.1 


0.56 ± 0.08 


1019.1 


0.59 ± 0.05 


1020.0 


0.51 ± 0.07 


1021.0 


0.74 ± 0.11 


1021.9 


0.81 lb 0.14 


1022.7 


0.55 ± 0.21 


1028.0 


0.53 ± 0.12 


1033.7 


0.97 ± 0.19 


1039.8 


0.96 ± 0.25 


1060.0 


1.08 ± 0.31 



Table 4: Measured cross section a{e~e'^ —>■ lut:^ tt^tt^j) (mean values between Cut I 
and Cut II, corrected by the factor 1.04, as discussed in the text). Only statistical errors 
are indicated. Systematic errors are estimated to be about 5%. 
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Figure 1: Theoretical predictions for cr(ee — > ujtt^ 7r°7r°7). Dashed hne - p-meson 
contribution only. Solid line: a) 0-meson contribution through — >■ pTT*^ — > 7r°7r°7 is 
added, b) — > ujt:^ transition due to p-oj-cj) mixing is also included in the 0-meson 
contribution, c) p'-meson contribution is added to the above ones. 



28 




20 40 60 80 100 cm 

Figure 2: Detector SND — view along the beam; 1 — beam pipe, 2 — drift chambers, 3 — 
scintillation counters, 4 — fiber lightguides, 5 — PMTs, 6 — Nal(Tl) counters, 7 — vacuum 
phototriodes, 8 — iron absorber, 9 — streamer tubes, 10 — 1 cm iron plates, 11 — scintillation 
counters, 12 — magnetic lenses, 13 — bending magnets 
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experiment 

Figure 3: Distributions M24 versus M35 plus M25 versus M34 for MC simulation and 
experiment. 
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Figure 4: versus xl-OnO-y distributions for MC simulation and experiment. 
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Figure 5: M^o^ (invariant mass of 7r'^7 nearest to uj) distribution for Cut I. 
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Figure 6: M^-o^ distributions : MC simulation (points with error bars) and experiment 
(histogram), a) for Cut I, and b) for Cut II. 
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Figure 7: Detected cross-section for the process e+e — > un^ — > 7r°7r°7 under cuts of 
group 1 and optimal fit. The fitted resonant background is also shown at bottom. 
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